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Abstract 

Excitonic spectra are calculated for free-standing, surface passivated InAs 
quantum dots using atomic pseudopotentials for the single-particle states and 
screened Coulomb interactions for the two-body terms. We present an anal- 
ysis of the single particle states involved in each excitation in terms of their 
angular momenta and Bloch-wave parentage. We find that (i) in agreement 
with other pseudopotential studies of CdSe and InP quantum dots, but in 
contrast to k.p calculations, dot states wavefunction exhibit strong odd-even 
angular momentum envelope function mixing (e.g. s with p) and large valence- 
conduction coupling, (ii) While the pseudopotential approach produced very 
good agreement with experiment for free-standing, colloidal CdSe and InP 
dots, and for self-assembled (GaAs-embedded) InAs dots, here the predicted 
spectrum does not agree well with the measured (ensemble average over dot 
sizes) spectra. (1) Our calculated excitonic gap is larger than the PL measure 
one, and (2) while the spacing between the lowest excitons is reproduced, the 
spacings between higher excitons is not fit well. Discrepancy (1) could result 
from surface states emission. As for (2), agreement is improved when account 
is taken of the finite size distribution in the experimental data, (iii) We find 
that the single particle gap scales as R~^-^^ (not i?^^), that the screened 
(unscreened) electron-hole Coulomb interaction scales as i?^^'^^ (R^^''^), and 
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that the eccitonic gap scla«s as i? These scahng laws are different from 

those expected from simple models. 
PACS:73.20.-r, 73.20.Dx, 85.30.Vw 
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I. INTRODUCTION 



Semiconductor quantum |T]] dots have recently generated considerable interest due to the 
dramatic differences between their electronic structure and that of the bulk material from 
which they are derived. These differences arise from the lower symmetry of the dot, quan- 
tum confinement level shifts and the enhancement of electron hole Coulomb and exchange 
interactions. 

While the Stranski-Krastanov growth technique produces semiconductor-embedded dots 
10-^ which are typically subject to only a small confining potential, colloidal chemistry 
techniques f^-Q produce dots whose surfaces are passivated by large organic molecules. 
This leads to much larger confining potentials (a few eV), and thus many more confined 
energy levels. These colloidal techniques have recently enabled the production of quantum 
dots with diameters from 10 to 60A , made from CdSe [Q, CdS 0, InP [T^ and InAs 



11|. The spectroscopy of these systems is richer than that of the SK dots, with several 



groups reporting data for up to 10 excited states in CdSe 0, InP [|T0| and InAs [|n[|. On 
the theoretical side, colloidal CdSe ||,|l3 and InP |[T0| , p!3| have recently been studied using 
both the standard 6x6 k.p method [§,[13] and more sophisticated pseudopotential |[T^ , p!3|JT^ 
techniques. 

In this paper we turn our attention to coUoidally grown InAs quantum dots. InAs is a 
challenging system to study because, (i) the small bulk band gap of 0.42 eV suggests strong 
valence-conduction band coupling. This is supported by recent studies [|ll|!0 showing that 
the standard 6x6 k.p formalism fails to describe the observed electronic transitions, while 
8x8 k.p produces a better fit to measured excitonic transitions, (ii) the small band gap 
combined with a large spin-orbit splitting (0.38 eV) suggests that spin orbit effects will play 
a significant role, (iii) the observed band gap extends to 1.6 eV, i.e. almost four times the 
bulk value, suggesting dramatic quantum confinement, (iv) the large confining potential 
leads to as many as eight clearly resolved excitations, for a range of dot sizes. 

We use a multi-band pseudopotential Hamiltonian to calculate both the single particle 
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states and electron-hole excitation energies for InAs quantum dots with a range of sizes. 
Our method is different from another recent pseudopotential calculation of Mizel and Cohen 



[ 151 in that although their formalism can be generalized to multiple bands, they use only 
a single-band approximation (like the EMA or the Truncated Crystal Method |T^) and do 
not include electron-hole effects. Both methods however, do incorporate non-parabolicity 
of the bands. The nature of the single particle states is analyzed in terms of the parent 
Bloch orbitals and envelope functions, which enables a detailed description of the origin of 
each calculated exciton. This analysis allows comparisons to be made with existing effective 
mass based predictions for the single particle states in these dots. In contrast to envelope 
function based effective mass calculations, our more general treatment shows that the single 
particle states include a significant amount of mixing of valence and conduction states, as 
well as a mixing of envelope functions with odd and even parity. We have calculated the 
single dot absorption spectra of a few dots. Comparing with the experimental results for 



ensembles of dots sizes, we do not find good agreement in contrast to the k.p method |Tl 
that gives better agreement. However, agreement with our results is improved when the 
finite size distribution of the experimental dot samples is taken into account. We discuss 
possible reasons for the conflicts between experiment and theory. 

II. PSEUDOPOTENTIAL CALCULATIONS OF THE ELECTRONIC 
STRUCTURE OF INAS QUANTUM DOTS 

Our pseudopotential calculations consist of two steps. First we calculate the single 
particle states for the quantum dot (Sections. [II A| and P B|) . Then we solve the "two- 



particle" problem by calculating the energy of the electron-hole excitations of the system 
(Section. 
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A. Single Particle Calculation 



We use a pseudopotential Hamiltonian to model the single-particle electronic structure 
of the system. 

H=-lv' + J2Mr-Ran)+vi'''^ . (1) 

The system's potential is constructed from a sum of screened atomic pseudopotentials, Va, 
where a represents In and As, and Rq,„ are the positions of the In and As atoms within the 
dot. The pseudopotentials, Va, have been fitted to the experimental band gaps and effective 
masses of bulk InAs. The experimental and fitted values are given in Table |. The fitted 
InAs bulk band structure is shown in Fig. [l|. We use the reciprocal space functional form of 
the pseudopotential 

q = aoa —2 , (2 

a2a&^°'^ — 1 

where aoo, ctia, a2a, «3a are adjustable parameters. The fitted values of ^ are given in 
Table The same pseudopotential form was recently used to study self-assembled GaAs- 
covered InAs pyramidal dots |T7|^0|. In Refs. |[r7|-pO| the pseudopotential contained an 
additional term to describe the effects of strain in the system. As the free-standing InAs 
dots studied here are strain free this term is not required here. We assume the bulk zinc- 
blende structure and the bulk In- As inter-atomic spacing (d=6.057A). We construct spherical 
InAs dots of radius i?, by selecting all the atoms that fall within a sphere of this radius. 
Any atoms from the surface of the dot which have three dangling bonds are removed. The 
result of adopting this atomistic description is to reduce the symmetry of the dots from 
the full spherical symmetry of continuum models to the lower T^ symmetry. The sizes and 



compositions of the four dots studied in this paper are listed in Table III 



To simulate the chemical passivation via ligands [0], we embed the InAs dots within an 
artificial barrier material, represented by an atomic pseudopotential, fitted to have a larger 
band gap than InAs, thus producing a type-I alignment between the dot and the barrier. The 
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barrier material is designed to have the InAs lattice constant so that no strain is introduced 
into the system. This embedding process is physically equivalent to the choice of a finite 
barrier surrounding the dots and is designed to refiect the fact that the dots are surrounded 



by organic molecules which themselves have large, but finite band gaps . 

We expand the single particle wavef unctions, ipi, in a plane wave basis whose cutoff must 
be identical to that used in the original pseudopotential generating procedure, 

Ecut 

^.(r) = ^ CG,e^«-^ . (3) 

G 

The matrix elements of the Hamiltonian in Eq.(0) in the basis of Eq.@ are calculated 
according to 

Hg,G' = 2^'^^G,G' + VlocaliG — G') + VnonlocaliG, G') . (4) 

The spin orbit interaction is represented by a non-local pseudopotential, VnoniocaiiG , G') , 
which is evaluated in real space using the linearly scaling small box method from Ref. ||22|| . 
This method applies the non-local pseudopotential to each atom in turn. For each atom, a 
new wavefunction ipbox{T^) = V^j(r) is defined within a small box around the atom, ipboxi^^) is 
periodically repeated and then fast Fourier transforms are used to generate ipboxiG). The 
non-local pseudopotential is then applied by 

UG)=J2ynonlocal{G,G') Aox{G') . (5) 
G' 

The InAs dots, surrounded by barrier material, form a supercell which is periodically 
repeated. Sufficient barrier atoms are used, to ensure that the interactions between an InAs 
dot and its periodic images are negligible. The total number of atoms (In, As and barrier) 
in the largest supercell is 25,000 atoms, which is too large for the Hamiltonian in Eq.(|l]) 
to be solved by conventional diagonalization methods. We thus use the "Folded Spectrum 
Method" (FSM) |2^J2^, in which one solves for the eigenstates of the equation 

[H - tref) ijji = {e - ereffi^i , (6) 
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where e^e/ is a reference energy. By placing e^e/ within the band gap of the dot, and close 
to the valence band maximum or conduction band minimum, one is then able to obtain 
the top few valence states or the bottom few conduction states respectively by minimizing 
< ip\{H—Eref)'^\ip >■ The use of FFTs results in an In scaling with the number of atoms, 
N, enabling calculation of systems containing 10^ atoms. It is easy to guess an approximate 
value of eref- We first determine the position of the bulk InAs band edges by using the same 
pseudopotential to calculate the bulk InAs band structure, e„fc. This is illustrated in Fig.|l|. 
As quantum confinement effects push electron(hole) levels up(down) in energy, placing eref 
within the bulk band gap ensures that it will also be within the band gap of the InAs dot. 
Note, that although Eq.(^ has twice as many eigensolutions as the standard Schrodinger 
equation, Hipi = ei ipi, we only search for solutions to Eq.(^ in a finite basis of A^ plane waves 
[Eq.(^)]. We therefore only have A^ solutions to Eq.(|^). By applying the H operator twice, 
it can be shown that each of these A^ solutions is a solution of Hipi = ipi and thus all our 
solutions to Eq.(||) are also solutions of the standard Schrodinger equation. Note also that 
our approach to the single-particle problem includes [Eq.(|T|)] shape effects, interface effects 
and spin-orbit coupling in the Hamiltonian [Eq.(|I|)]. In the effective mass approximation, 
on the other hand, these effects are included perturbatively [|rT|. Also note that we neglect 
self-consistency with respect to the bulk. As the dots contain several hundred atoms we 
expect the potential of the dot interior to be bulk-hke. The pseudopotential used here was 
fitted to reproduce these bulk properties without requiring a self-consistent calculation. 



B. Methods for analyzing the single particle wavefunctions 

Having calculated the single particle states of an InAs quantum dot, we are interested 
in analyzing the nature of each of these states in terms of the bulk Bloch wavefunctions 
{0n/c} from which they are derived, the symmetry of their envelope functions and their 
defining quantum numbers. These quantities will enable us to classify each of the electron- 
hole excitations calculated in the following section in terms of the detailed nature of their 
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initial valence state and the final conduction state. They also enable us to make contact 
with alternative electronic structure theories, such as the k.p technique, which typically 
classifies states in terms of their envelope function and total angular momentum. Note, the 
expressions in Eqs.(|^-(|1^) below are used only to analyse the solutions of Eq.(^) and do 
not affect the results. 

We first deconstruct each single particle dot wavefunction, ipi{r), in terms of the bulk 
Bloch states at the T point, 0„r = Mnr(i")e*''''" in the manner described in Ref. [jl2|. We choose 



this basis to (i) enable the calculation of the different angular momentum character of the 
states, and (ii) to enable comparison with conventional envelope function based methods. 
These zone center Bloch states form a complete set, and therefore one can expand the single 
particle dot wavefunction according to 



^ik.r 



^n,r(r) (7) 



U^) = E 

n 

which can be rewritten as 

^,(r) = 5:«„,r(r)/«(r) , (8) 

n 

where u„^r(r) is the bulk Bloch wavefunction for the nth. band at the T point, and /^*^(r) is 
a corresponding envelope function. Note that the "single band approximation" corresponds 
to retaining a single n value in Eq.(§), while the 6x6 k.p corresponds to n = bulk VBM 
(sixfold degenerate, including spin). The present "multiband" method corresponds to a 
large number of bands in Eq.(|^). We further decompose each envelope function, /^*^(r), into 
spherical harmonics, FL^m(6',0), according to, 

/«(r) = E<U(|r|)>^L,a^,</>) , (9) 

L,m 

and then define the weight from each angular momentum component, w^\, as 



n,L 



E / \rRlmir)\'dr , (10) 



where R is the radius of the quantum dot. The quantity, w^ ^.) tells us how much of bulk 
band n and angular momentum L is contained in the dot wavefunction ipi. 
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An additional property that is useful for classifying the single particle states is the total 
angular momentum, F, of dot state i. This is calculated for each of the single particle dot 
states according to 

i^(^+i) = (Ewi^^i(^«+^/)'iEwi:A • (11) 

\ ra ml 

The angular momentum operator F„ = ir x V + s, (where s is the Dirac spin matrix), acts 
only on the Bloch orbitals, M„(r), and = x V acts only on the envelope function, 

C. Two particle calculation and the absorption spectrum 

Calculation of an optical absorption spectrum based on the single particle states of a 
quantum dot, requires a solution for the "two-body" electron-hole problem. More complete 
solutions of the electron-hole problem via the configuration interaction require knowledge of 



the Coulomb, exchange and correlation energy associated with the electron-hole pair 
In this work, we choose to ignore the fine structure effects resulting from the exchange and 
correlation energy, reducing the two-body problem simply to calculating the electron-hole 
Coulomb energy. We then define the exciton energy in terms of the initial single particle 
valence state energy, e^, the final conduction state energy, e^, and the electron-hole Coulomb 
energy, Jj^, as 

^exciton — i^j ~ ^i) ~ Jij ■ (12) 

The electron-hole Coulomb energy, Jij, between each of the possible electron(i) - hole(j) 
pairs is calculated by a direct integration of the single particle wavefunctions 

where ipi{r) and are our calculated electron and hole wavefunctions and e(ri — r2,R) 

is a distance dependent screening function, connected to the dielectric function, e(r,r',i?), 

by, 

9 



and n is a size scaling exponent discussed below. We follow Ref. [0 in writing the Fourier 
transform of e^^(ri — R) as the sum of electronic and ionic terms, 



e~\k;R)=e;i\k;R) + AerJ^{k;R) 



(15) 



The electronic term is defined in terms of the Thomas Fermi model of Resta 26 



Q\k;R) 



e + q^ sm{kR^)/{et'kR, 



(16) 



Note, the R dependence of e^i^{k; R) comes from e^*. We define q in Eq.(^6|) as 



(17) 



where Hq is the electron density, and R^o is the solution to 

sinh(gi?oo) 



qRc 



(18) 



The ionic term is taken from Ref. IE? 



1 



1 



1/2 



+ 



1/2 



(19) 



Here ph,e = i'^'^'h e'^Lo/^)"^^^; ^Lo is the longitudinal optical-phonon frequency (238.6 cm^^) 
and mj^g are the electron and hole effective masses (0.023 and 0.4 mo). To obtain the value 
of eo°*, we again assume that the dot interior is bulk-like so eg"* — e^* can be approximated 
from the bulk, 



^dot ^dot ^bulk ^bulk 



-0 



-oo ~ "=0 



Ae: 



bulk 



(20) 



We use the bulk value of 2.9 for Ae^"^'^. Finally, we use a generalization |^ of the Penn 
model PP[ to obtain ej^*. 



et\R) = 1 + (e 



bulk 
oo 



(« + A)2 



(21) 
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where A = E2 - E^"^^ = 4.28 eV. The electronic and ionic components of the screening 



gap 



function are plotted in Figs. 0(a) and (b). One sees that the electronic contribution dom- 
inates for all but the smallest A;-vectors. Figure ^ shows the screened dielectric function 
in real space, e{R,r), as a function of electron-hole separation r, for dots with radii, R, 
from Table At very small separations, e{R,r) approaches 1.0, which corresponds to no 
screening at zero separation between the electron and hole. When the electron and hole are 
separated by r = 6 — 8A, e(r) approaches the particular value of e^*(i?) which is different 
for each of the dots : The resulting values of e^*(i?) are 7.19, 7.91, 8.41 and 8.92 for dot 
diameters 23.9, 30.3, 36.6 and 42. 2A. The ionic screening, Ae^on has a very long tail, such 
that e(r) approaches eo°* for r > lOOA. The R-dependence of e(r, R) can be clearly seen in 
Fig. ^. In all cases the value is smaller than the bulk InAs value of eg"''^ = 15.15. 

More simplified calculations of J^- assume a universal value for all i and j, based on an 



s-like envelope function and an infinite potential barrier . These approximations lead to 
the often adopted formula 

jEMA ~3.572 , 

Our pseudopotential calculated values for the screened electron-hole Coulomb energy be- 
tween an electron in the lowest energy conduction state and a hole in the highest energy 
valence state, Jn, are plotted along with those predicted by Eq. (p2D in Fig. ^. The magni- 
tude of the pseudopotential values is larger than that predicted by Eq.(]22|) for all sizes of 
dot. A fit of the dependence of the screened Coulomb energies, J, on the dot radius, R, 
to J = A/R"-, yields a size exponent, n, of 1.79, for the pseudopotential results, compared 
to the 1.0 scaling predicted by the Eg. p^) . The difference in magnitude between Eq.(p2D 
and pseudopotential calculated Coulomb energies of Eq. (|13]) has two sources: Firstly, we use 
microscopic (rather than envelope) functions that are not required to vanish exactly on the 
surface of the dot as EMA functions are required to do when subjected to an infinite poten- 
tial barrier. This wavefunction effect causes the unscreened Coulomb energy eJ to scale as 
with m < 1, while in the EMA m = 1. This effect was discussed by Franceschetti et. 
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al. who found m ~ 0.7 for GaAs. Secondly, the description of screening in Eqs. (|T^) to 
leads to a reduced screening at small R with respect the bulk screening, e'^''^. This acts 
to increase our values of Jij and its R dependence with respect to those in Eg. (p2D . This 
effect is oppostie and stronger than the wavefunction effect. The net effect is J ~ i?"^ ''^ 
compared with J ~ R^^ in Eg. (P^D . The present calculated values for the size scaling of the 
unscreened electron-hole Coulomb interaction, eJ, and the screened Coulomb interaction, J 



for InAs dots are compared with calculations for Si |^ , InP ||38| and CdSe [|T^ dots as well 



as the predictions from Eg. (^) in Table It shows that for the scaling of eJ has m < 1 
while the scaling of J has n > 1. 

The final stage in the calculation of our theoretical absorption spectra, is to calculate 
the dipole matrix transition element for the transition from valence state i to conduction 

state j, Mi^j, 

M,^^{R) = mV\^P,)f . (23) 
The optical absorption spectra, I{E,R), can then be written as. 



/(E,i?) = ^^M,_,(i?)exp 



exciton 

a 



(24) 



where E^^l^^^^ is the excitation energy [Eg.(|T^)] and a represents an experimental linewidth 
broadening which we set to 25 meV. 

Our theory of Eg.(|^) corresponds to the spectra of a single dot with a well defined 
radius, i?, and shape. Current samples exhibit a size distribution about an average size, i?, 
producing an ensemble absorption spectra 

1{E, <R>) = j P{R)I{E, R)dR (25) 

where P{R) is the size distribution function. A major contribution in the interpretation of 
experimental data is that if the size distribution is too broad, a given transition {ij) for size 
Ri could overlap with another transition (i'j') for size R2- 
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III. RESULTS 



A. Single Particle States 

Using the Hamiltonian described in Section jlT^ , the 48 highest energy hole states and 



the 24 lowest energy electron states were calculated for the InAs dots listed in Table pi . 
Figure |^ shows plots of the electron and hole wavefunctions squared, for the dot with a 
diameter of 42. 2A, in a plane through the center of the dots. The plots show that as a 
result of the large band offsets between the vacuum and the InAs dot, the wavefunctions 
are strongly localized within the InAs dot. Approximately 90.5% and 95.9% of the ground 
state electron and hole wave function squared respectively, are localized within the dot. 

The single-particle energies of the highest energy hole state and the lowest energy electron 
state are plotted as a function of diameter in Fig. |^. For comparison, the same states obtained 
from a recent Wannier function based, single-band pseudopotential calculation of Mizel and 
Cohen are also shown. Figure ^ shows that in both calculations, as the size of the InAs 
dot increases, the quantum confinement decreases, such that the VBM and CBM approach 
the bulk VBM and CBM values. We have fitted the size dependence of the dot VBM and 
CBM according to 



dot 
vbm 


jpbulk 

vbm 


+ 


a 


dot 
cbm 


Tpbulk 

cbm 


+ 


b 



(26) 

where R is the radius of the InAs dot. This fitting reveals that the size exponents, n„^c, of 
the hole and electron quantum confinements are significantly different to the effective mass 
exponents: 



0.95 Current Pseudopotential 

1.23 Single band pf ■ (27) 

2.0 EMA 



1.21 Current Pseudopotential 
nv = <, 1.59 Single band [15] ric = 

2.0 EMA 

Figure |^ also shows that while the CBM states calculated with the multi-band pseudopoten- 
tial and the single-band method are almost identical in energy, the VBM states calculated 
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using the single-band method are higher in energy than those calculated using the pseu- 
dopotential. This difference can be attributed to the fact that the real dot valence band has 
strong interband coupling, so the single band approximation used by Mizel and Cohen 



is more severe for the valence band. For the CBM states, the single band approximation 
is sufficient because the energy spacing between bands is large and there is therefore little 
interband coupling. 

Single Particle Gaps: In Fig. 0(a) we plot the single particle band gaps obtained from our 
multi-band pseudopotential calculations. For comparison we also show the single particle 



gaps obtained from the single band calculations [15| of Mizel and Cohen and the 8 band k.p 



calculations from Ref ||32|. Figure ^(a) shows that the single band calculations |jT5| tend to 



underestimate the single particle gap with respect to the more converged calculations, as a 
result of the underestimate of the valence band quantum confinement in single band models. 
The 8 band k.p and current multi-band pseudopotential results show good agreement with 
each other. In Fig. ^(b) we compare the multi-band pseudopotential calculations with a 
recent PL experiment, by extracting the single particle gaps from experimental results from 
Ref. 10 by subtracting an approximate value for the electron-hole Coulomb energy calculated 
from Ecoui = 3.572/2e-R (Ref. [§D|)- Here e is the static dielectric constant of bulk InAs and 
R is the dot radius in atomic units. We also show the results of the single-particle gap 
extracted from a new STM tunnelling experiment on the same colloidal InAs dots. The 



multi-band pseudopotential calculations predict larger gaps than those measured by PL. 
This can have three reasons (i) the observed emission could involve surface states below the 
CBM and above the VBM and would therefore have a lower energy than the calculated band 
to band values, (ii) The calculation assumes a single dot, while in the experiment there is 
an ensemble of dots. Thus, the observed emission could be of low energy since the largest 
dots in the ensemble, having the lowest CBM and highest VBM energies emit. In addition 
(3) measured sizes suffer from the notorious difficulty of measuring size accurately. If the 
measured TEM size is underestimated, (especially for the samll dots) this could explain 
the discrepancy. The new tunelUing STM experiments on similar InAs dots, Banin et. al. 
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measure a single particle gap that is higher than that extracted from the optical experiment 
and is in much better agreement with the calculations. 

We fit the single particle (sp) band gap to the following power law, 

E J* = E,^;^''^ + AD" , (28) 

taking the bulk gap as 0.42 eV. Our pseudopotential results yield a size exponent Ugp = 
— 1.01. This is considerably smaller than the effective mass predicted value of -2. In the 
previous section we have seen that the Coulomb interaction scales as R~^'^, while here 
we see that the single particle gap scales as Thus the relative importance of Coulomb 

effects increases for small dots. This is the opposite of the effect predicted by single-band 
effective mass predictions of J{R) ~ R^^ and Egap ~ R~^. The scaling, n, of the single 
particle band gap is compared with similar calculations for InP, CdSe and Si dots in Table 0. 

Excitonic gaps: In Fig. |^(c) we plot the actual PL results along with our multi-band pseu- 
dopotential excitonic results calculated from Eq.(|12[) including the electron- hole Coulomb en- 
ergy screened according to Eg. (p!3|) . The single-particle pseudopotential results from Fig. |^(a) 
are shown for comparison. The size exponent for the excitonic gap is Uexciton = —0.90, smaller 
than that for the single particle gap (n^p = —1.01). Using the Coulomb energy from Eq.(|l3D 
rather than Eq.(^) improves the agreement with experiment in Fig. ^(b) by a factor of 3. 
However, the fit is still not particularly good. The reasons for this were discussed above. 

Confinement energies: We have also calculated the ratio of the conduction band quantum 
confinement, AEcb{R), to the valence band quantum confinement, AEysiR), defined as 

AT? { D\ ^1^"* ^bulk 

AEvb{R) = (^ttL - (^Ibm ■ (29) 

Define ^(-R) as 

Our multi-band calculations for InAs dots find a value of ^{R) that decreases with R. The 
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best linear fit (constant ^(-R)) to the data yields a ratio, ^(-R), of 3.13 for the multiband 
calculations. For Si dots, C,{R) was recently measured ||3^ to have a value of 0.5. 

Within the single-band effective mass approximation, assuming an infinite potential bar- 
rier at the surface, the quantum confinement of electrons and holes can be written as 
Aeg./j = ^^7r^/2m* ^i?^, where m*^ is the effective mass of electrons and holes at the F 
point and R is the size of the quantum dot. The ratio of electron confinement to hole 
confinement 

^{R) = — = — = const. 31) 

is therefore given by the ratio of the electron and hole effective masses which is approximately 
15.4. This is significantly different to our multiband calculated value of 3.13. This difference 
can be attributed to the fact that the pseudopotential calculated hole states are derived 
form a mixture of the bulk heavy and light hole states. 

B. The nature of the single particle wavefunctions 

The results of the single particle wavefunction analysis from Eqs.(|8|)-(pilD is given for the 
42. 2A diameter dot in Table |V|. Details of the analysis for the other dot sizes are given at 



Ref. 1^ . Only the states which contribute to the excitonic peaks discussed in Section |1I C 
are listed in the table. For each state, the fraction of the wavefunction derived from the 
bulk F-like split off, heavy (Fg^) and light hole (Fy^,) and lowest conduction band (Fgc) states 
is given (see Fig. ||). For each of these bulk states, the fraction of the total wavefunction 
derived from envelope functions with s, p and d symmetry is also given. The significant 
contributions to each state are marked in bold. Contributions less than 0.01 have been set 
to zero. For example, the highest energy hole state has a fraction of 0.0 + 0.02 + 0.02 = 0.04 
derived the bulk split-off band, 0.50 + 0.29 + 0.09 = 0.89 derived from the bulk heavy and 
light hole bands and 0.0 + 0.02 + 0.0 = 0.02 derived from the lowest bulk conduction band. 
The remaining 0.05 fraction is derived from bulk bands further from the band gap, and from 
higher angular momentum envelope functions. 

16 



Analysis of the results in Table and those for the other size dots reveals several inter- 
esting properties of the single particle wavef unctions: 

(1) The origin of the lowest lying electron states in the dot follows qualitatively the 
predictions of single band effective mass theory. For example, the lowest electron state of 
the 42.2A diameter dot (see Table [V]) is 69% derived from the bulk conduction band edge 
state (Fee) with an s-like envelope function. The next two highest electron levels are 64% 
derived from the same bulk Bloch state, but with a p-like envelope functions. 

(2) As a result of the small band gap of bulk InAs, there is a strong coupling between 
the electron and hole states. Approximately 20% of the weight of the lowest energy elec- 
tron states in the dot is derived from the split-off, heavy and light hole (Fg^ + Tj^)) bulk 
states. This valence-conduction mixing explains why the 6x6 k.p method, which ignores 
such coupling, fails to describe these states in InAs dots. 

(3) The highest energy hole states in the dot have significant weight from both the s and 
p envelope functions from both heavy and light hole bulk states. They therefore cannot be 
described, even qualitatively, by a single band model. The VBM of the 42. 2A diameter dot 
(see Table |V]) has 50% s and 29% p character. Similarly, the lowest electron state in the 
dot has 27% non-s character, originating from valence bands. Such an s-p mixing is largely 



absent in current theoretical descriptions of InAs dots via the k.p method |Tl| . 

(4) The order of the electron and hole states changes with size. This reflects different size 
scaling of the quantum confinement for different states. For example in the 42. 2A diameter 
dot (see Table 0) the lowest conduction state is s-like (69%), then states 2,3 and 4 are p-like 
(64%), then states 5,6,7 and 8 (not shown in the table) are a mixture of d and g. In the 
30. 3A diameter dot, the states with 42% p character have moved above the states identified 
by their significant d and g character making them states 6,7 and 8 as opposed to 2,3 and 
4. Similarly, in the valence band the state with 37% s-like contribution from the split off 
band and 13% and 12% from p and d envelopes of heavy and light hole states is the 5*^ and 
6*'' state from the top of the valence band in the 30. 3A diameter dot but moves to the ll*'^ 
and 12*^ state from the top of the valence band in the 42. 2A diameter dot. 
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C. Calculated Single-dot absorption spectra 



Theoretical single-dot absorption spectra were calculated for each of the dots listed in 
Table |T| using the method described in section |II C| . These spectra are shown in Figs. |^(a)- 



(d). The identities of each of the major peaks in the spectra were determined by examining 
the nature of the initial and final single particle states contributing to each peak. The 
following criteria were used to establish these identities. 

1. The ratio of the contributions to the single particle states from the heavy hole, light 
hole, split off and conduction states. 

2. The contribution of each angular momentum component from the heavy hole, light 
hole, split off and conduction states. 

3. The total angular momentum, F, of each single particle state (see section [11 B| ). 

4. The strength of the dipole transition probability matrix element for each peak. 

Table |V| shows the identities of the peaks for the 42. 2A diameter dot. Details of the analysis 
for the other dot sizes are given at Ref. They show that there are two separate manifolds 



of transitions: peaks (a) to (e) which correspond to transitions from the dot hole states to 
the lowest s-like electron state and peaks (f) to (k) which correspond to transitions from the 
dot hole states to the next highest p-like electron state. The analysis of the peaks shows: 

(a) "peak a" represents the fundamental band gap transition. The initial valence state 
associated with this peak is a doubly degenerate hole state with s and p character derived 
from both heavy and light hole bulk states. This initial state has a total angular momentum, 
F, that ranges from 1.75 to 2.01. The final conduction state is the singly degenerate, lowest 
energy electron state, with an s-like envelope function and a total angular momentum that 
ranges from 1.11 to 1.27. In the approximate k.p language, this transition is closest to the 
5*3/2 to Si/2 transition. 

(b) "peak b" corresponds to a transition with very weak intensity. For the largest dot 
with 42. 2A diameter, this state has merged with peak (a) and cannot be resolved. The initial 



valence state associated with this peak is a doubly degenerate hole state with p character 
derived from both heavy and hght hole bulk states. It is the second highest energy hole 
state. This initial state has a total angular momentum that ranges from 1.92 to 2.51. The 
final conduction state is the same as in peak (a). In the approximate k.p language this 
transition is closest to the P5/2 to S1/2 transition. 

(c) "peak c" has a strong intensity. The initial valence state associated with this peak 
is a doubly degenerate hole state with a mix of s, p and d character. Its most significant 
contribution is from the split-off state but it contains some heavy and light hole character. 
This initial state has a total angular momentum that ranges from 1.99 to 2.48. The final 
conduction state is the same as in peak (a). In the approximate k.p language this transition 
is closest to the P5/2 to S1/2 transition. 

(d) "peak d" is a transition with very weak intensity. For the largest dot with 42. 2A di- 
ameter this transition merges with "peak c" and is no longer distinguishable. The initial 
valence state associated with this peak is a doubly degenerate hole states with a mix of 
mostly s and d character. It is derived from the bulk heavy and light hole states. This ini- 
tial state has a total angular momentum that ranges from 2.48 to 3.25. The final conduction 
state is the same as in peak (a). In the approximate k.p language this transition is closest 
to the 5*5/2 to 5*1/2 transition. 

(e) "peak e" has a similar origin to "peak d" and also has a very weak intensity. It 
also merges with "peak c" in the largest dot with 42. 2A diameter. The initial valence state 
associated with this peak is a singly degenerate hole state with a mix of s, d and some g 
character. It is derived from both the bulk spin-orbit, heavy and light hole states. This initial 
state has a total angular momentum that ranges from 3.03 to 3.75. The final conduction 
state is the same as in peak (a). In the approximate k.p language this transitions is closest 
to the 5*7/2 to 5i/2 transition. 

(f) "peak f" corresponds to a transition with very weak intensity that is only observed 
in the two smaller dots. It has the same initial state as "peak a" , but the final state is the 
next highest conduction state. This conduction state is a triply degenerate state, with a 
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p-like envelope function and a total angular momentum that ranges from 2.15 to 2.19. In 
the approximate k.p language, this transitions is closest to the 6*3/2 to P3/2 transition. 

(g) "peak g" corresponds to a transition from the same initial state as "peak b" , to the 
p-like conduction state. This is a very intense peak that is observed in all the sizes of dot. 
In the approximate k.p language, this transitions is closest to the P5/2 to P3/2 transition. 

(h) "peak h" is also a very intense peak that is observed in all the sizes of dot. The initial 
valence states are a combination of the same initial state as "peak c" and also a doubly 
degenerate hole state, very close in energy to this state with mostly s and / character, 
derived from heavy and light hole states, with a total angular momentum ranging from 
2.66 to 2.72. The final conduction state is the same as in peak (f). In the approximate 
k.p language, this transition is closest to the P5/2 to P3/2 transition. It is possible that 
the proximity of peaks (g) and (h) combined with their strong intensity would not allow 
them to be distinguished in a photoluminescence experiment. This is discussed further in 
section [1IIIl| . 

(i) "peak i" is has a weaker intensity, but is observed in all the sizes of dot. The initial 
valence state is the same as "peak d". The final conduction state is the same as in peak 
(f). In the approximate k.p language, this transition is closest to the S5/2 to P3/2 transition. 
This transition is weakly allowed. 

(j) "peak j" corresponds to a transition with strong intensity. The initial valence state 
associated with this peak is a doubly degenerate hole state containing mostly a mix of p, d 
and / character. Its contains significant contributions from both the bulk split-off state and 
the heavy and light hole states. This initial state has a total angular momentum that ranges 
from 3.03 to 3.41. The final conduction state is the same as in peak (f). n the approximate 
k.p language, this transition is closest to the 5*7/2 to P3/2 transition. 

(k) "peak k" corresponds to a transition with weak intensity. The initial valence state 
associated with this peak is a doubly degenerate hole state containing mostly p character. 
It is mostly derived from the bulk heavy and light hole states. This initial state has a total 
angular momentum that ranges from 3.19 to 3.24. The final conduction state is the same 
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as in peak (f). In the approximate k.p language, this transition is closest to the P7/2 to P3/2 
transition. 

Having discussed the identities of the peaks in each dot size, we next wish to see how to 
connect peaks with the same identities in different sizes of dot. This is not always possible, 
as different sizes of dot might have some peaks that are fundamentally new, or two peaks 
that have merged together. We have labelled in Fig. || and in Table the peaks that 
originate from similar excitons in the different sizes of dot by the same letters (a) to (k). 

D. Comparison of calculated single-dot spectra with the experimental absorption 

spectra 

The positions of the above peaks are plotted as a function of diameter for each dot in 
Fig. ^ and with respect to the band gap of each dot in Fig. |10[ As some of the peaks merge 



with each other in the larger dots, not all the peaks are marked for all the sizes of dot. 
The predictions of Figs. |^ and |10| pertain to hypothetical samples each containing dots of 
a single size and spherical shape. The actual synthesized samples |jTl| contain a significant. 



but unknown, distribution of dot sizes and shapes. Nevertheless, the experimental results 
from Ref. [|rT| are included in Figs. ^ and |TU[ It should be noted that where two peaks 
have effectively merged together, or where the weight of a particular peak is too small to be 
detected (e.g. peak "b"), these points are not plotted in Figs. |^ and |IU[ A comparison of 
the pseudopotential results in Fig.^ with the results from Ref. [jll| is given in Table |VI1| , 
which shows that 

1. Peak (a), by definition, corresponds to the experimentally measured band gap. El. 

2. We ascribe two weak sets of peaks(c and d) as originating from the weakly observed 
E2 experimental peak. 

3. For the E3 peak we calculate a single peak (e). 

4. For the weakly observed E4 peak we find no calculated counterpart. 
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5. The strongest two calculated peaks (g and h) fall on either side of the E5 peak and it is 
possible that the strength of these excitations could prevent them from being isolated 
in the experiment. They are also merged by a finite size distribution (see section ^TT^). 

6. The final two calculated peaks (j and k) correspond to the experimental E6 and E7 
data. 

For each of the peaks, the calculated scaling of the exciton energy spacings with dot size 
(or band gap) shows reasonable agreement with the experimental results. However, Fig. ^ 
shows that the calculated values of the absolute exciton peaks appears to exhibit a different 



size dependence to that observed in Ref. |TT[. This lack of agreement can be attributed in 



part to the finite size and shape distribution present in the experimental samples. This is 
discussed in the following section. 



E. Ensemble absorption spectra 

The ideal comparison between theory and experiment is between the calculated (Fig. ^ 
and measured single-dot spectra. However, no such single-dot measurements currently exist 
for InAs quantum dots. Our predicted single-dot spectra for different sizes (Fig. |]) suggest 
that the interpretation of an ensemble spectra could differ qualitatively from the interpre- 
tation of a single dot spectra. This calls for measurement of the single dot spectra. For 
example, peak (b) in the dot with a diameter of 30.3 A coincides with peak (e) in the dot 
with a diameter of 36.9 A and with peak (g) in the dot with a diameter of 42. 2A. Thus, if 
the experimentally accessible samples represent a broad size distribution it is impossible to 
ascribe consistently experimental peaks to unique calculated single-dot peaks. 

We try to quantify the effect of a finite size distribution in the experimental samples by 
using our single dot spectra in combination with Eq.(p4D to calculate ensemble absorption 
spectra. These are not directly comparable with the size selected PLE results from Ref. 
0, but are designed to provide a general indication of the effects of size distribution on 
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the absorption spectra of an ensemble of dots. We neglect shape distribution effects, since 
they are not quantified experimentally. Transmission electron microscopy studies of III-V 
semiconductor quantum dots show that there are two factors producing an ensemble 
of different dot volumes. Firstly, in any sample there is a finite range of dot diameters. 
Secondly, the dots are ellipsoidal in shape, with a range of ratio of major to minor axes. 
We have therefore chosen to model the distribution of sizes, P{R), in Eq. (p4|) by a simple 
Gaussian, whose width, a^, builds in the size distribution, 

P(i?) = ^^e-(^-^«)'/2<xi _ (32) 
^/2^TaJi 

In Fig. |ll|(a) we plot ensemble absorption spectra calculated from Eq.(|^) for quantum dots 
with a mean diameter of 23.9 A and standard deviations, a^, of 0, 5 and 10% of the mean 
size. The function I{E, R) in Eq.(]2^) was obtained by fitting the size dependence of each of 
the peaks, i, in Fig. ^ to Ei{R) = Ef + ai?^ and then summing the contributions from all 
the peaks so that 

IiE,R)= J2 a.e-(^-^»(^»'/2.g ^ (33) 

peaks,i 

where is the relative intensity of peak i, and uo is the intrinsic line width of the peaks, set 
to 5 meV. Figures 0(b) to (d) show similar ensemble absorption spectra for dots with size 
distributions that have mean diameters of 30.3, 36.9 and 42.2 A and standard deviations 



of 5%. Figures |TT|(a) to (d) show that as the width of the size distribution increases, the 
peaks become "smeared out". With a standard deviation of 5% it is still possible to resolve 
peaks, (a), (c), (g) and (h). However, peak (j) is just a small shoulder on peak (h). With a 
standard deviation of 10% peaks (g) and (h) have effectively merged together. The merging 
effect becomes more severe as the mean size of the dots increases as the initial spacing of 
the peaks is smaller for larger dots. 

Figure ^ shows how the centers of the broadened peaks in each of the 5% ensemble 
spectra from Fig. |Tl] vary with size. Figure |12] is the size ensemble equivalent of the single 
dot results plotted in Fig. |TU|. As in Fig. 0, the experimental results of Banin et. al. are 
plotted for comparison. Comparison of Figs. and |12|, reveals that 
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1. The fundamental transition (a) is still clearly resolvable in the ensemble spectra and 
corresponds to experimental peak Ei. 

2. The first excited state (c) is still clearly resolvable in all sizes of dot. This peak 
probably still corresponds to peak E2. 

3. Peak (d) is only resolvable in the two smallest dots, where it is close to the experimental 
E3 peak. 

4. For the larger two dots, peaks (g) and (h) merge to form one large peak. The position 
of this merged peak is close to the experimental E5 peak. For the two smaller dots, 
the peak splits into two peaks (g) and (h) with different size scaling behaviour. 

5. Peak (j) is only resolvable from peaks (g) and (h) in the two smaller dots, where it 
could correspond to either the experimental Eg or E7 peaks. 

6. The weaker peaks (b), (e) and (k) are not individually resolvable for any size of dot in 
the ensemble spectra. 

IV. CONCLUSIONS 

We have performed pscudopotential calculations for the electronic structure of both the 
ground and excited states of free standing InAs quantum dots for a range of experimentally 
realistic sizes. Using calculated electron-hole Coulomb energies and dipole matrix transi- 
tion probabilities we have constructed single-dot absorption spectra for 4 different sizes of 
quantum dot. These spectra exhibit a series of clearly resolvable exciton peaks. The size 
dependence of the spacing between the exciton peaks in these single dot spectra shows par- 
tial agreement with those found in recent experiments. By fitting the size dependence of 
each exciton peak and postulating a Gaussian distribution of dot sizes, we have calculated 
ensemble absorption spectra. The size scaling of the peaks in these ensemble spectra shows 
a better agreement with the experimental data. 
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We also have analyzed the single particle parentage of each excitonic peak. We find 
that (i) as a result of the small band gap of InAs, there is significant valence-conduction 
band mixing in the quantum dot states, (ii) the removal of spherical symmetry of these dots 
produces odd-even mixing in these states. 

Our predicted single-dot excitonic spectra (Fig. ^) await experimental testing. Our 



predicted ensemble spectra (Figs. |Tl|) are not in as good agreement with experiment as our 
results for InP [0] and CdSe [|l2l dots. Sample quality, including shape distributions could 
be a factor in this relative lack of agreement. 
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TABLES 

TABLE L Experimental and fitted bulk InAs pseudopotential properties. Aso is the spin-orbit 

splitting and m^^^, "ipis^ hh anFig. 0(a)d ^^fi5„ ih ^he effective masses of electrons, heavy and 
light holes. 

Property Expt. Value |37| Pseudopotential 

Tic - ri5^ (eV) 0.42 0.41 

Xic - Ti5^ (eV) 2.33 2.27 

Lie - (eV) 1.71 1.61 

Ao (eV) 0.38 0.36 

ml 0.029 0.028 

A Ic 

mK,.,hhi'^00] 0.43 0.41 

mf^^^^^JlOO] 0.038 0.039 
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TABLE II. Screened atomic pseudopotential parameters for InAs and the barrier potentials. 



The a parameters refer to Eq.(|2|). 

ao ai 02 a-i 

In 118.781 1.783 3.382 0.393 

Barrier Cation (cond) 107.755 1.915 3.460 0.414 

Barrier Cation (val) 333.070 0.120 3.126 0.521 

As 65.943 2.664 1.684 0.637 

Barrier Anion (cond) 19.892 2.097 1.182 0.243 

Barrier Anion (val) 49.614 2.737 1.523 0.574 
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TABLE III. InAs quantum dot sizes and compositions 



Dot Number 12 3 4 

No. In atoms 140 276 456 736 

No. As atoms 141 249 459 683 

Diameter (A) 23.9 30.3 36.6 42.2 
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TABLE IV. A comparison of the size scaling exponent, n, in of the unscreened elec- 
tron-hole Coulomb energy ,eJ, the screened electron-hole Coulomb energy, J, the single particle 
band gap, EgPp and the excitonic band gap E^^^^""" [Eq.(|l2|)] for InAs, InP, CdSe and Si quantum 



dots. 



gap 



Property 



InAs 



InP HI 



CdSe 121 



Si HI 



EMA 



eJ 
J 



gap 

rpexciton 
gap 



0.70 
1.79 
1.01 
0.90 



1.35 
1.16 
1.09 



0.86 
1.18 
1.64 
1.93 



0.82 
1.50 
1.20 
1.18 



1 
1 
2 
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TABLE V. Contributions (as fractions of unity) from the split-off (SO), heavy hole (hh), light 
hole (Ih) and conduction band (CBM) bulk states to the single particle states of the InAs dot with 
42. 2A diameter. The dot states are numbered (first column) from the band edge. For electron 
states 1 is the lowest in energy, 8 the highest. For hole states, 1 is the highest in energy, 24 the 
lowest. The contribution from each state is determined using Eq.(|8|). For each state the bulk SO, 
hh+lh and CBM contributions are decomposed into their angular momentum components using 
Eq.([lO|). Only the states contributing to the major excitonic peaks in Fig. ^ are shown. All entries 
less than 0.01 are set to 0.0. The main contributions are marked in bold. Totals include envelope 
functions with angular momenta from to 6. 



Level 


Energy (eV) 


SO Contribution 


hh + Ih Contribution 


CBM Contribution 


s 


p 


d 


Total 


s 


V 


d 


Total 


s 


p 


d 


Total 


Conduction States 


























1 




0. 


0.07 


0. 


0.07 


0. 


0.20 


0. 


0.20 


0.69 


0. 


0. 


0.69 


2 


€c + 0.360 


0.06 


0. 


0.01 


0.07 


0.01 


0. 


0.19 


0.20 


0. 


0.64 


0.01 


0.65 


3,4 


ec + 0.361 


0. 


0. 


0.07 


0.07 


0.09 


0. 


0.10 


0.19 


0. 


0.64 


0.01 


0.65 


Valence States 


























1,2 




0. 


0.02 


0.02 


0.04 


0.50 


0.29 


0.09 


0.89 


0. 


0.02 


0. 


0.02 


3,4 


6^,-0.014 


0. 


0.01 


0.01 


0.03 


0.23 


0.55 


0.07 


0.88 


0. 


0.01 


0. 


0.01 


7,8 


- 0.098 


0. 


0.02 


0.02 


0.08 


0.08 


0.38 


0.18 


0.78 


0. 


0. 


0.02 


0.02 


11,12 


6^,-0.123 


0.28 


0. 


0. 


0.31 


0. 


0.26 


0.19 


0.56 


0. 


0. 


0.01 


0.01 


13 


- 0.140 


0. 


0.21 


0.06 


0.28 


0. 


0.12 


0.21 


0.60 


0. 


0. 


0. 


0. 


15,16 


-0.175 


0. 


0.06 


0.01 


0.14 


0.02 


0.26 


0.05 


0.73 


0. 


0. 


0. 


0. 


18 


- 0.189 


0. 


0.04 


0.13 


0.27 


0.07 


0.16 


0.04 


0.69 


0. 


0. 


0. 


0. 


23,24 


- 0.223 


0. 


0.05 


0.01 


0.14 


0.01 


0.33 


0.12 


0.69 


0. 


0. 


0.02 


0.02 



33 



TABLE VI. Analysis of the significant absorption peaks for a dot with 42. 2A diameter. The 
peak letters refer to the labelling of the peaks in Figs. |l^ and 0. Peak energies are calculated 



according to Eq. ( 12 ) . The valence and conduction indices refer to the number of the state from the 
band edge (i.e. VBM and CBM have index 1 etc). The envelope function information summarizes 
the results from Table ^ Note, the notation "s", "p" etc refers to the orbital content of the 
wavefunctions, not the specific number of nodes. All energies are in eV. 



Peak Peak Energy 



Initial valence state 



Final conduction state 







Index 


Energy 


Envelope 


F 


Index 


Energy 


Envelope 


F 


(a) 


1.209 


1,2 


-5.8260 


s + p 


1.77 


1 


-4.516 


s 


1.11 


(c) 


1.335 


11,12 


-5.9491 


p + d 


2.48 


1 


-4.516 


s 


1.11 


(g) 


1.588 


3,4 


-5.8401 


P 


1.90 


2,3,4 


-4.156 


P 


1.92 


(h) 


1.675 


7,8 


-5.9240 


s+P + f 


2.66 


2,3,4 


-4.156 


P 


1.92 


(i) 


1.717 


13 


-5.9660 


p + d + f 


2.93 


2,3,4 


-4.156 


P 


1.92 


0) 


1.760 


15,16 


-6.0087 


P + f 


3.07 


2,3,4 


-4.156 


P 


1.92 


(k) 


1.800 


23,24 


-6.0491 


P 


3.19 


2,3,4 


-4.156 


P 


1.92 
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TABLE VII. Comparison of the current multi-band pseudopotential assignment and the k.p 
assignment from Ref. [^] of the experimentahy observed optical transitions in Ref. |11|. The initial 
states and final states of each transition are denoted by nQp, where n is the principal quantum 
number, Q the lowest angular momentum component of the wavefunction and F total angular 
momentum. Note, it is only possible to determine the envelope function angular momentum and 
total angular momentum for the multi-band pseudopotential wavef unctions. 



Peak 



Pseudopotential assignment 



Closest Expt. Peak 



from Ref. |11| 



k.p assignment 



from Ref. |11| 



(a) 


5*3/2 


Sl/2 


El 


25*3/2 - 


15*1/2 


(b) 


P5/2 


Sl/2 


Not observed 


No prediction 


(c) 


P5/2 


Sl/2 


E2/E3 


No prediction 


(d) 




Sl/2 


E2 


No prediction 


(e) 


S7/2 


Sl/2 


E3 


25*3/2 - 


15*1/2 


(g) and (h) 


P5/2 


^P3/2 


E5 


1^^3/2 - 


- 1^3/2 










1^^3/2 - 


lPl/2 


(j) 


S7/2 


P3/2 


Eg 


25*1/2 - 


15*1/2 


(k) 








1^^1/2 - 


- 1^3/2 


P7/2 


P3/2 


E7 










1^^1/2 - 


- 1^/2 
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FIGURES 

FIG. 1. In As bulk band structure from the L to F to X calculated with the empirical pseu- 
dopotential defined in Eq.(|2|). 

FIG. 2. The electronic and ionic contributions to the dielectric screening function. Plotted in 
reciprocal space. 

FIG. 3. The screened dielectric function. Plotted in real space for each of the 4 InAs dots given 



in Table HI. 



FIG. 4. Electron- hole Coulomb energies calculated using both pseudopotential [Eq.(|l3|)] and 
effective mass [Eq.(|2^)] expressions for all 4 sizes of InAs quantum dot given in Table p^ . 



FIG. 5. Highest energy hole and lowest energy electron wavefunctions squared in a plane 
through the center of the dot for a free standing spherical InAs quantum dot with diameter of 
42.2 A. 

FIG. 6. Size dependence of the lowest energy hole state and highest energy electron states, 
calculated using the present multi-band pseudopotential and single band techniques. The circles 
and triangles mark the present multi-band pseudopotential and single band (Wannier function) 
results. The solid lines are the best fit to a scaling of Eq + AR"" (see Eq.(p6|)). 

FIG. 7. (a) Calculated single particle band gaps calculated using present multi-band pseu- 
dopotential (circles), 8x8 k.p (solid line) and single band techniques (squares), (b) Measured single 
particle band gaps using Photoluminescence (solid line) and STM experiments (triangles), (c) 
Present multi-band pseudopotential results (circles) including the electron-hole Coulomb energy 
calculated from Eq.(|l^) vs PL data. 

FIG. 8. Calculated absorption spectra for single InAs quantum dots with diameters of (a) 
23.9, (b) 30.3, (c) 36.6 and (d) 42. 2A. The labels (a)-(k) refer to the classification of the states in 
Table Ivl. 
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FIG. 9. Pseudopotential calculated (circles) absolute positions of absorption peaks vs 
for InAs quantum dots with diameters 23.9, 30.3, 36.6 and 42. 2A. The pseudopotential results are 
connected using the analysis from section III Q . Experimental results from Ref. are represented 
by the black lines. 

FIG. 10. Pseudopotential calculated (circles) positions of absorption peaks minus the band gap 
vs band gap for InAs quantum dots with diameters 23.9, 30.3, 36.6 and 42. 2A. The pseudopo- 



tential results are connected using the analysis from section IIIC. The labels (a)-(k) refer to the 



classification of the states in Table VI. Experimental results from Ref. |11] are represented by the 
black symbols and lines. 

FIG. 11. Calculated absorption spectra for ensembles of InAs quantum dots with mean diam- 
eters of (a) 23.9, (b) 30.3, (c) 36.6 and (d) 42. 2A. (a) shows ensemble spectra for size distributions 
with 0, 5 and 10% standard deviations, (b), (c) and (d) show only the spectra associated with a 



size distribution. The labels (a)-(k) refer to the classification of the states in Table Vl. 



FIG. 12. Pseudopotential calculated (circles) positions of absorption peaks minus the band gap 
vs band gap for InAs quantum dot ensembles with mean diameters of 23.9, 30.3, 36.6 and 42. 2A and 
standard deviations of 5%. The pseudopotential results are connected using the analysis from 
section III C. The labels (a)-(k) refer to the classification of the states in Table VI. Experimental 



results from Ref. |11] are represented by the black lines. 
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